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ABSTRACT 

We describe the infrared properties of sources detected over ~36 deg 2 of sky in the GAM A 15-hr 
equatorial field, using data from both the Herschel Astrophysical Terahertz Large- Area Survey (H- 
ATLAS) and Wide-field Infrared Survey (WISE). With 5ct point-source depths of 34 and 0.048 mJy 
at 250 /im and 3.4 /im, respectively, we are able to identify 50.6% of the H- ATLAS sources in the 
WISE survey, corresponding to a surface density of ~ 630 deg -2 . Approximately two-thirds of these 
sources have measured spectroscopic or optical/near-IR photometric redshifts of z < 1. For sources 
with spectroscopic redshifts at z < 0.3, we find a linear correlation between the infrared luminosity 
at 3.4 /im and that at 250 /im, with ±50% scatter over ~ 1.5 orders of magnitude in luminosity, 
~10 9 - 10 10 - 5 L Q . By contrast, the matched sources without previously measured redshifts (r > 20.5) 
have 250-350 ^m flux density ratios that suggest either high-redshift galaxies (z > 1.5) or optically 
faint low-redshift galaxies with unusually low temperatures (T < 20). Their small 3.4-250 /im flux 
ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low 
rcdshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large 
AGN fraction (~30%) in a 12 /im flux-limited subsample of H- ATLAS sources, also consistent with 
there being a significant population of high-redshift sources in the no-redshift sample. 
Subject headings: surveys - cosmology: observations — galaxies: high-redshift - infrared radiation - 
galaxies: statistics - galaxies: general 
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1. INTRODUCTION 

The advent of submillimeter astronomy has opened a 
new window into the Universe, allowing us to probe dusty 
star-forming galaxies at high redshift, as well as cold 
dust in nearby galaxies. Along with the recent success 
of the balloon-borne B LAST experiment (Pasc ale et al.l 
l2008HEales et al.ll2009l), the launch of the Herschel Space 
Observator$\ ()Pilbratt et al.l l2010f ) allows us to probe 
wavelengths from 55 to 672 /iin to a spatial resolution 
of < 10". Its largest open-time key project, the Herschel 
Astrophysical Terahertz Large- Area Survey (H- ATLAS, 
lEales et ai1l2010al) . will observe — 550 deg 2 of sky and 
detect more than 300,000 galaxies. 

Based upon a preliminary cross-identification with 
the Galaxy and Mass Assembly survey (GAMA, 
Driver et al.l[20ll an d Sloan Digital Sky Survey (SDSS, 
Abazaiian et al.l l2009). -50% of the sources detected in 
the H- ATLAS survey are at z < 1 (|Smith et al.ll2011l) . In 
addition, an analysis of the FIR colors of the remaining 
H-ATLAS sources with flux densities above 35 mJy at 
350 /im and > 3cr detections at 25 /zm and 500 (tm sug- 
gests an average redshift of z — 2 (|Amblard et al.ll2010l ; 
ILapi et a.l.lf201lTl . 

Another method of constraining the redshift distri- 
bution of FIR-selected sources is to look for counter- 
parts in the near and mid infrared (NIR and MIR). 
Fleuren et al. (in prep) have performed source match- 
ing to survey data from the VISTA Kilo-degree Infrared 
Galaxy survey (VIKING, Sutherland et al., in prep), but 
these data only extend to — 2 /im. Starting in December 
2009, the Wide-fie ld Infrared Survey ExploreiQ ( WISE, 
Wrig hTet al.l 120 lOh began its mission to observe the en- 
tire sky in four bands, ranging from 3.4 to 22 /im, at 
< 12" resolution. Following the first public data release 
in April 2011, — 24,000 deg of IR images and source cat- 
alogs became available to the public, including 36 deg 2 
of sky in the equatorial plane covered by the H-ATLAS 
survey. 

The purpose of this letter is to describe the in- 
frared properties of identified WISE counterparts to LI- 
ATLAS sources within the GAMA 15-hr field (G15). 
Throughout we will assume a concordance cosmology 
with ff n = 71 km s" 1 Mpc" 1 , Q m = 0.27, and Q A = 0.73 
(jSpergel et al.ll2007h . In a subsequent paper, we will re- 
port WISE cross-identifications over the entire H-ATLAS 
area and perform spectral energy distribution (SED) fits 
to the matched sources. 

2. DATA AND METHODOLOGY 

Herschel observations of the G15 field include imag- 
ing data at 250, 350, and 500 /im from the SPIRE in- 
strument, wjth^esEectiyebeam FWHM of 18"1, 25"2, 
and 36"6 (| Griffin et al.ll2010f ). as well as 100 and 160 /zm 
imag ing from the PACS instrument (iPoglitsc h et al.1 
2010). The field subtends approximately 12° in right 
ascension and 3° in declination and is centered on the 
equatorial plane at a = 14 h 30 m. The H-ATLAS G15 
source catalogue (Dunne et al., in prep) is constructed 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA. 

2 For a description of the WISE m ission, see 
http : / / wise2 .ipac . caltech.edu / docs/release / prelim / expsup / 



using the same method as the Sci ence Demonstration 
Phase catalogue (|Rigbv et al.l 120111 ). It contains 27,481 
sources detected at > 5cr in any of the three SPIRE bands 
and reaches point-source depths of 34, 40, and 44 mJy 
at 250, 350, and 500 /xm, respectively. 

The WISE first public data release contains four-band 
coverage of the entire G15 field to 5a point-source depths 
of 0.048, 0.10, 0.73, and 5.9 mJy at 3.4, 4.6, 12, and 
22 /im. The angular resolutio n in these bands is 6"1, 6'.'4, 
6"5, and 1270, respectively (| Wright et al.l[20Tot ). There 
are ~ 240,000 WISE sources within G15 detected at > 7a 
in at least one of the four WISE bands, corresponding to 
— 0.14 WISE sources per 10"-radius aperture. In the 
WISE preliminary release catalog, close pairs of sources 
are not deblended for separations < 9", so when a back- 
ground or foreground source appears near an H-ATLAS 
source position, it will often be blended with the true 
counterpart. 

In addition to the infrared data from WISE and H- 
ATLAS, the G15 region has spe ctroscopic redshift s from 
the GAMA survey (r < 19.4. iDriver et al.ll20lil ), and 
photometric redshifts obtained using optical/NIR pho- 
tometry from SDSS, VIKING, and the UK Infrared 
Deep Sky Survey L arge Area Survey (UKIDSS-LAS, 
lLawrence et al J 120071). These p hotom etric redshifts are 
derived following iSmith et al.1 ()201lD and have typical 
redshift uncertainties of a z /z — 0.15. Of the 27,481 H- 
ATLAS sources in the region, 16% have reliable spectro- 
scopic redshifts and 39% have photometric redshifts. 

As a result of the high space density of 3.4 /im sources, 
a naive matching to the H-ATLAS source positions will 
result in a non-negligible number of misidentifications. 
A significant fraction of these misidentifications will be 
foreground stars that are blue in all WISE bands and 
undetectable in H-ATLAS, so we can reduce the con- 
tamination rate of our matched catalog by considering 
only WISE sources with [3.4] - [4.6] > 0. The magnitude 
distribution of all background/foreground WISE sources 
is such that — 20% are bluer than this limit, as compared 
to < 1% of sources within 10" of an H-ATLAS position. 

To estimate the fraction of H-ATLAS sources with 
detectable counterparts in the remainder of the WISE 
source catalog, we use the method of Fleuren et al., in 
prep, where the detection rate is given by, 

Qo = i-§, (i) 

where S is the fraction _of unmatched positions in the 
H-ATLAS catalog and B is the fraction of unmatched 
random positions. Matching all WISE sources within 
10" (within which we expect > 99% of the true matches 
to lie) and using 10 5 random positions within the G15 
field, we find Q = 0.632 ± 0.004. This is larger than 
the value found in the SDSS galaxy catalog (Qo = 0.583, 
ISmith et al.ll20TlT) . but smaller than in the VIKING Ks 
band (Qo = 0.75, Fleuren et al., in prep). By contrast, 
we find Q = 0.012 ± 0.002 within the subset of WISE 
sources with [3.4] — [4.6] < 0, suggesting that this color 
cut was effective in removing stars. 

In Fig. [U we use Eq.[T]to estimate the IR detection rate 
as a function of FIR flux density. For sources brighter 
than 100 mJy at 250 /im, we find a WISE 3.4 /im coun- 
terpart > 95% of the time. The majority of these sources 
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Fig. 1. — NIR and MIR detection rates plotted as a function 
of limiting FIR flux density at 250 fim (solid), 350 (im (dashed), 
and 500 /im (dotted). Detection rates are shown for 3.4 /mi (upper 
left panel), 4.6 (im (upper right), 12 /im (lower left), and 22 ', fim 
(lower right) sources, where the corresponding 5cr detection limits 
are 0.06, 0.9, and 3.6 mjy, respectively. Note that the large sky 
density of 3.4 /mi sources can lead to ambiguity in the identification 
of some sources, so identification rates are typically ~ 10-20% lower 
than the corresponding detection rates. 



(78%) have measured spectroscopic or photometric red- 
shifts < 0.5, suggesting that the low-redshift mode of 
star-forming galaxies dominates sub-mm sources above 
this flux density. Below this FIR flux density, the NIR 
and MIR detection rates drop rapidly as the high-redshift 
mode accounts for an increasing fraction of the sub-mm 
sources (see Fig. [5]). 

Approximately 1% of the ob j ects in our 250 fim- 
selected sample (jNegrello et al.l I2007T) . including ~ 
50% of sources brighter than 100 mj y at 500 /im 
(jNegrello et alj2010tlHopwood et al1l2011l ). are expected 
to be strongly-lensed z > 1 galaxies. Of the 46 such 
sources in the G15 region, ~20% are undetected at 12 
and 22 /im, most likely because the 500 /im bandpass 
shifts higher in the blackbody curve at high redshift, 
while the rest- frame MIR flux d ensity declines blueward 
of 22 /im (e.g. JRieke eTal][2009l ). 

The identification of individual NIR and MIR coun- 
terparts is more subtle, as we want to eliminate as 
many of the false matches as possible. Here, we use 
the lik elihood ratio technique of Sutherla nd fc Saunderd 
(HHl), which, was i mplem ented, for the H-ATLAS sur- 
vey in ISmith et all (|20lU hereafter Sll). The likeli- 
hood that a given WISE source is a counterpart to an 
H-ATLAS source is a function of the radial probability 
distribution, /(r), the 3.4/im magnitude distribution of 
non-counterparts, n(m), and the 3.4 /im magnitude dis- 
tribution of detected H-ATLAS counterparts, q(m): 



q(m)f(r) 
n(m) 



where /(r) is given by 



i 
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Fig. 2. — Distribution of the 3.4— 250 (im flux density ratio as a 
function of 250 /im flux density for H-ATLAS sources with WISE 
3.4/mi IDs. Colored points have spectroscopic redshifts from the 
GAMA survey and sources without measured redshifts are indi- 
cated by black points. Dashed lines indicate the approximate lim- 
its of the H-ATLAS and WISE surveys. The majority of sources 
with spectroscopic redshifts (r < 19.4) are well detected by WISE, 
while the no-redshift sources have small flux density ratios that are 
inconsistent with "normal" low-redshift galaxies. 

We measure n(m) directly from the WISE prelim- 
inary release catalog, while we determine q(m) by 
first measuring the magnitude distribution of 3.4/im 
sources within 10" of the H-ATLAS source positions 
and then subtracting the magnitude distribution of 
background/foreground sources normalized by the area 
within this radius. 

Following Sll, we assume the astrometric uncertainty 
to be circularly symmetric and to have a dependence on 
the 250 fim signal-to- noise ratio (SNR) given by, 



FWHM 
Tpos = °- 655 ^NR-' 



(4) 



where FWHM= 18.1" is the Herschel beam width. While 
the WISE random astrometric uncertainties are negli- 
gible (~0"3), the WISE positions in the preliminary 
release catalog are known to suffer from highly non- 
Gaussian systematic offsets that can be as large as 1", so 
we apply a lower limit of cr pos > 1.5". 

Once we have computed the likelihood, Lj, for a pos- 
sible counterpart j, we determine the reliability by sum- 
ming over all possible counterparts, 



L, 



+ (i-Qo) 



(5) 



We consider any 3.4 /xm source with R > 0.8 to be a 
reliable counterpart. After performing this procedure 
separately for WISE sources with star- and galaxy-like 
[3.4] — [4.6] colors, we obtain a combined total of 13,898 
WISE counterparts, or 50.6 ±0.4% of the H-ATLAS cat- 
alog. The expected number of false matches can be ob- 
tained with 
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We estimate that there are 369 false identifications in 
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Fig. 3.— FIR flux density ratio for H- ATLAS sources with WISE 
identifications. Colored grid cells indicate the median redshift of 
a set of modified blackbody models with temperatures uniformly 
distributed between 10 and 60 K and emissivity parameters be- 
tween and 2. Black points indicate sources with 5cr detections 
in all three SPIRE bands, where open squares have spectroscopic 
redshifts, z < 0.2, and solid triangles have no spectroscopic or 
photometric redshift measurement. Dashed lines indicate the me- 
dian flux density ratios of a stack of sources lacking redshifts. The 
majority of the no-redshift sources appear to be at z ~ 2. Inset: 
Combined distribution of spectroscopic and photometric redshifts 
for H- ATLAS sources with detections at 3.4 fim (solid histogram), 
12 /im (dashed), and 22 (im (hatched). There are 73 sources at 
z > 1 not shown. 

our sample, corresponding to a contamination rate of 
2.7%. The majority of our IDs have galaxy-like colors, 
but we were able to identify 14 objects with star-like 
colors ([3.4] — [4.6] < 0) that weren't already identified 
as galaxies with SDSS matching. This corresponds to an 
upper limit of 0.05% for the fraction of H- ATLAS sources 
identified as stars and detected in WISE. 

Comparing the individual WISE identifications to 
SDSS cross-identifications, we find that 84.7% of the 
10, 709 G15 H- ATLAS sources with R > 0.8 SDSS identi- 
fications also had R > 0.8 WISE identifications. The me- 
dian separation between the WISE and SDSS sources is 
tf.'6, consistent with the astrometric uncertainties in the 
current version of the WISE catalog. In addition, ~ 6% 
of these sourcess have separations > 3", consistent with 
the expected contamination rates of the WISE (~ 3%) 
and SDSS (~ 5%) matched samples. 

3. DEMOGRAPHICS 

The redshift distribution of 250 /im-selected sources 
is likely bimodal, with populations of both moderately- 
star-forming galaxi es at z < 1 and high-redshift starburs t 
galaxies at z ~ 2 (iSmith et al.l l2011t lLapi et alJfeOllD . 
Because the latter are typically faint in the observed- 
frame optical and UV, the vast majority of objects with 
measured spectroscopic or photometric redshifts are in 
the low-redshift mode. Of the H-ATLAS sources with 
spectroscopic and photometric redshifts (see the inset of 
Fig.©, we successfully identify 85.6% in WISE. 

Although most of the WISE identifications with known 
redshifts are at z < 0.8 (>98%), such a sample is 
biased toward low redshift because redshift measure- 



ments require detectable flux in the UV, optical, or near- 
infrared. Approximately one-third (30.4%) of the iden- 
tified sources have reliable spectroscopic redshift and 
two-thirds (66.5%) have photometric redshifts from op- 
tical/NIR data. For the remaining sources, we can esti- 
mate redshifts using their FIR flux density ratios (e.g., 
lAmblard etldl 12010: iSchulz et all 12010 ). which give an 
indication of the location of the peak of the FIR dust 
emission. In Fig. [3l we compare the SPIRE colors of our 
matched catalog to a suite of 10 6 modified black body 
spectral energy distributions (|Amblard et al.ll2010f ) . We 
indicate on the plot the median flux density ratios of the 
4, 677 H-ATLAS sources with WISE identifications but 
no redshifts ("no-redshift" sample, dashed lines), as well 
as the individual flux density ratios of 311 such sources 
with > 3cr detections in all three SPIRE bands. 

The no-redshift sample has median FIR colors sug- 
gesting z ~ 2, while the subset with three-band SPIRE 
detections may have even higher redshifts, perhaps as 
large as z ~ 3.5. The redshifts derived using the 
lAmblard et al.l (|2010[) technique are highly uncertain, as 
redshift is degenerate with dust temperature and it is 
possible that these galaxies have typical dust temper- 
atures diff erent than 35 K (whi ch is typical of high-z 
SMGs, e.g. lChapman et al"1 f2005). However, in order for 
the no-redshift sample to be at z ~ 0, they would need 
to have very cold dust temperatures (T ~ 12 K) in ad- 
dition to being very faint in the optical (r > 20.5). Al- 
though temperatures as low as 10 K have been seen in 
H-A TLAS galaxies a t z < 1, T ~ 25 - 30 K is more typ- 
ical ()Dve et al.ll2010f ). If the no-redshift sample has me- 
dian temperatures typical of these low-redshift H-ATLAS 
galaxies, then our redshift estimate drops to z ~ 1.5. 

Further evidence that the no-redshift sources are at 
z ^> can be found in the 3.4-250 /im ratios (see Fig. [2]). 
H-ATLAS sources with GAMA redshifts have NIR-FIR 
flux density ratios that decrease from ~ 0.009 at z ~ to 
^0.006 at z ~ 0.35. This trend, which continues toward 
higher redshift, is due primarily to the larger fc-corr ection 
at 3.5 /im than at 250 /im (e.g.. iRieke et al] 12009). The 
median flux density ratio for the no-redshift sources, by 
comparison, is 0.0023. Although low-redshift galaxies 
will occasionally exhibit such small NIR-FIR flux density 
ratios - Arp 220, for example, has /3.4//250 = 0.0019 - 
such objects are rare, actively star-forming, and unlikely 
to have dust temperatures as small as ~ 15 K. 

We can shed even further light on the properties of 
the no-redshift sample by examining their position in 
WISE color space (see Fig. @|. The technique is sim- 
ilar to those developed for a set of IRAC filters at 
comp arable wavelengths (|Saiina et al]|2005l : iStern et all 
[2005h . It was demonstrated that the diagram can be 
divided into stars and early-type galaxies (lower left cor- 
ner, [4.6] — [12] < 1), star- forming galaxies (lower right 
corner), and active galactic nuclei (AGNs), where AGNs 
are selected largely based on their power-law emission 
in the NIR and M IR. This was adapted for WISE by 
Uarrett et all (|2011ft . 

[4.6] - [12] > 2.2 
[4.6] - [12] < 4.2 
[3.4] - [4.6] < 1.7 (7) 
[3.4] - [4.6] > 0.1Q4.6] - [12]) + 0.38. 
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Fig. 4.— IR color distribution of ATLAS- WISE matches with 
three-band WISE photometry, including sources with z < 0.5 
(color coded by redshift) and no-redshift sources (black points). 
The diagram can be divided into stars (lower left corner), star- 
forming galaxies (lower right corner), and AGNs (dashed region). 
The color distribution of all three-band WISE detections (including 
sources unmatched to H- ATLAS) is indicated by contours. Sources 
without measured redshifts tend to be redder in [3.4] — [4.6] and 
have a larger AGN fraction than sources with spectroscopic or pho- 
tometric redshift measurements, z < 0.5. 

Of the 4,959 H-ATLAS sources with WISE identifica- 
tions, z < 0.5, and detections in 3.4, 4.6, and 12 /mi, the 
majority are star-forming galaxies, with < 0.2% main 
sequence stars (which is consistent with the SD SS/H- 
ATLAS matching done bv lThompson et alJfeOlOt ). and a 
small AGN fraction (0.057), 

For the 9% of no-redshift sources with three-band 
WISE photometry, we find a much larger AGN frac- 
tion (0.30), but are subject to a selection bias toward 
AGNs due to their addition al flux from warm dust at 
12 /im llHainline et al.l [201 If) . Furthermore, emission 
from polycyclic aromatic hydrocarbons in high-redshift 
star-forming galaxies can occasionally mimic AGNs in 
their IRAC c olors, leading to contamination in our AGN 
samples (e.g . lLacv et al.ll2004t iDonlev et al.l l2012D. For 
comparison, iCoppin et al.l (|2010l ) performed SEP fits on 
a set of blank-field SMG samples and estimated that 
~ 15% of SMGs are dominated by an AGN in the mid- 
infrared. 

4. FIR-NIR PROPERTIES 

Emission from dust in the FIR is often used as an in- 
dicator of the to tal star formation rate in galaxies (e.g., 
IKenmcuttlll998D . under the assumption that most of the 
dust heating is provided by young stars and that the 
star-forming regions are optically thick. The 250 /im 
bandpass samples the cooler dust emission, which can 
arise from bot h star-forming clouds and diffuse regio ns in 
the ISM fe.g. JEales et~aT1l2010btlDunne et alJl201lh . By 
contrast, 3.4 /im emission is dominated by stellar sources 
and is a probe of the stellar m ass for galaxies wi thout a 
strong AGN contribution (e.g.. iStern et al.ll2005T) . 

In Fig. [SI we plot 3.4 /zm luminosity as a function of 
250 /xm luminosity for H-ATLAS sources with spectro- 
scopic redshifts, 0.05 < z < 0.3. We fc-correct 3.4 /im 
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Fig. 5. — Luminosity at 3.4 /im as a function of luminosity at 
250 fj,m for H-ATLAS sources with spectroscopic redshifts, 0.05 < 
z < 0.3. Both luminosities are rest-frame fc-corrected. Small 
squares are colored according to their redshift and red stars in- 
dicate AGN, selected with the criteria given in Eq. [7] Dashed lines 
indicate the approximate selection limits for sources at z = 0.05 
and 0.2. The majority of the AGN lie above the relationship, sug- 
gesting that L3.4 contains a non-negligible contribution from hot 
dust emitting near the central black hole. 

using a power law interpolation of the [3.4] — [4.6] color. 
For the 250 /im flux densities, the majority of the galax- 
ies in our sample do not have sufficient FIR flux den- 
sity to obtain reliable dust temperatures, as is needed 
for a proper k correction. Instead , we use the me dian 
dust temperature of 26 K found by I Dye et al.l (|2010T ) us- 
ing PACS and SPIRE data from the H-ATLAS SDP of 
1,346 0.1 < z < 0.5 sources. Although the r < 19.4 limit 
of the GAMA spectroscopic survey restricts the range 
of luminosities that can appear on the plot, we estimate 
that GAMA misses only ~ 5% of matched z < 0.3 galax- 
ies based upon the deeper (r < 20.5) subsample of H- 
ATLAS galaxies with photometric redshifts. 

We find that the majority of the IR-color-selected AGN 
(red points) lie above the relationship for star-forming 
galaxies, as -L3.4 will contain a non-negligible contribu- 
tion from hot dust emitting near the central black hole 
in addition to emission from low-mass stars. However, if 
we exclude AGNs, we find an approximately linear cor- 
relation, with a best-fit power law index, a = 0.98±g g|. 
The quoted systematic uncertainties on the power law 
index were determined by allowing for a range of possi- 
ble dust temperatures, 18 < T < 34 K, when performing 
250 /im fc-corrections. The intrinsic scatter about this 
relationship is 0.18 ± 0.01 dex, or ~ 50%. 

The existence of a linear correlation between L250 and 
L3.4 across one and a half decades in luminosity suggests 
a close relationship between the cold dust probed by the 
FIR and the stellar mass probed by the NIR. Previous 
indications with Herschel have shown the 250 //m lumi- 
nosity density to tightly correlate with both the 24 /im 
luminosity density, a star for mation rate indicat or, and 
the total infrared luminosity (lElbaz et al.l |2010|) . This 
fact, coupled with a linear correlation between stellar 
mass and star form ation rate for the g e neral star-forming 
galaxy population (jDaddi et al J 120071: lElbaz et al.l 120071: 
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iNoeske et al.ll2007l Donoso et al., in prep), suggests that 
the majority of the low-redshift H-ATLAS sources are 
actually "normal" star-forming galaxies. 

This publication makes use of data products from the 
Wide-field Infrared Survey Explorer, which is a joint 
project of the University of California, Los Angeles, 
and the Jet Propulsion Laboratory /California Institute 
of Technology, funded by the National Aeronautics and 



Space Administration. 

Herschel is an ESA space observatory with science in- 
struments provided by European-led Principal Investiga- 
tor consortia with sig- nificant participation from NASA. 
U.S. participants in Herschel ATLAS acknowledge sup- 
port provided by NASA through a contract issued from 
JPL. 

In addition, we thank Dan Stern and Roberto Assef 
for helpful discussions about the quasar selection. 
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